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SUMMARY 


The, results of calculations of the viscous and pressure drags 
of some two-dimensional supersonic airfoils at zero lift are presented. 
The resuJLts indicate that inclusion of viscous drag alters ‘many 
previous results regarding the desirability of certain airfoil shapes 
for securing low drags at supersonic speeds . • At certain Reynolds and 
Mach nvuiihers, for instance, a clrculcir-arc airfoil may theoretically 
have less drag than the previously advocated symmetrical wedge -shape 
profile; althou^ under different conditions, the circu3ar-arc 
airfoil may have the higher drag. 

Drag calculations for 6 -percent -thick symmetrical circular-arc 
and double -wedge airfoils are presented for Mach numbers of 1.35 

and '1.6 and Reynolds numbers from 10^ to 10®, Unseparated flows are 
considered and approximate corrections for boundary -layer and shock - 
wave interaction are applied only to the momentum thiclmess at the 
trailing-edge shock wave. The iiieory of viscous supersonic flows 
will have to be extended before a more exact aneilysis of the drag is 
possible. 

INTRODUCTION 


Recent experiments indicate that airfoil shapes heretofore 
considered good for sxipersohic speeds may in fact be inferior to 
profiles haying higher pressure drags. In order to understand better 
the behavior of airfoils at supersonic speeds, it is desirable to 
eliminate and explain apparent contradictions betwen experiment and 
theory. The drag of thin airfoils laay be considered as the sum of the 
pressure and viscous drags., Althoxigh a great deal of work has been 
done pn ^ie. calculation of the pressure drag, little theoretical work 
has been attempted on the calculation of the viscous drag. The purpose 
of this paper is to consider the separate effects of viscous and 
pressure drag on 1±ie total drag for wo thin airfoil sections. 
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Hie pressure drag as determined"By"’tEe approxinjate linearized 
relations is a minimum for- a; sya^trioa^ \re^e-sha^ airfoil of a 
given thickness ratio i-riLth maximum thickness located at the 
50 "percent chord (reference 1). The; mo3?e -exacjj methods of refer" 
ences 2 and 3 show that other locations of the maximum thickness for 
the symmetrical i^edge "shape profile have. Ipiver pressure drag. On i 
these flat"side profiles, hoi'rever, the velocity gradient is zero over 
the flat leading and trailing sections' of the airfoils and, hence, 
"boundary "layer transition and viscous dra,g may he somewhat similar to 
that of a flat plate. 'Viscous* drags lower than those usually obtained 
for a flat, plate can he obtained by using a falling "pressure gradient 
to increase the extent of the laminar flow. It appears possible then 
that a ciurved airfoil, such as a double circular arc, which has a 
favorable pressure gradient may also have /more laminar flow than a 
flat "Side airfoil and therefore less viscous drag over a certain range 
of Mach and Eeynolds numbers. Although the theoretical pressvre drag 
of the circular "arc profile is higher, than that of the symmetrical 
wedge shapes, the total drag (viscous and pressure) of the circular 
arc might possibly be the loi-7er for certain conditions. 

In order to demonstrate the effects of viscous and pressure ■ drag 
on the total drag, calculations were made for two^ 6 "percent "thick 
airfoil shapes " a double circular arc and a symmetrical wedge with 
maximum thickness at the 50"P®^cent chord. The calculations were 
made for Mach numbers of 1.35 ^d 1.6 and covered a range of Eeynolds 
numbers from 10^ to 10°. 

The present paper is intended to serve only as a preliminary 
study of the total drag of a two "dimensional airfoil and neglects 
many factors. A more complete analysis would reqtiire consideration 
of the effects of separation, interaction between the boundary layer ' 
and the shock waves, and angle of attack. Existing theories need 
considerable development before all these factors can be included. 


DISCUSSION 


The drag of a t\ro -dimensional airfoil at zero lift at supersonic 
speeds is assumed to be the sum of the pressure drag and the drag due 
to viscosity. The pressure distribution. is caiculated by assuming 
the absence of the boundary layer, and the . shock drag is then readily 
determined from the pressure distribution* The .boundary -layer 
momentum thickness corresponding. ..to., this, pressure distribution is 
calculated, and the viscous drag is then determined from the momentum 
thickness at the trailing edge of the airfoil. 
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03ie method of reference 4 was used to determine the shock drag 
as well as the local values of Mach numher, velocity^ density, and 
other conditions along the airfoil. 03iis method requires the 
existence of attached shock waves and is iiierefore restricted to 
sharp*nose airfoils. In the absence of a "boundary layer and flow 
separation the calculations are accurate for the wedge airfoil and 
are a close approximation for the circular-arc airfoil. In a more 
complete analysis the pressure distribution should be adjusted for 
changes in the flow pattern caused by bovindary -layer thickness and- 
by .sudden changes in slope of the boundary -layer displacement 
thickness duo to transition or separation. 

The calculated pressure distribution is vised to coE^iute the 
boundary -layer momentum thickness along the airfoil by the method 
of reference 5 . Reference 5 assumes the following; The skin- 
friction coefficient is independent of !4ach nvmiber'and pressure 
gradient j a fixed velocity profile independent of pressure gradient 
may be used; the Prandtl number is 1 ; and no heat conduction occurs. 
Both laminar and turbulent boundary layers may be. computed approxi- 
mately by this method with the use of the. appropriate constants 
given in the reference. For the present calculations, transition 
from laminar . to turbulent flow was considered to occur suddenly. 
Therefore, in order to con^iute the momentum thickness along the surface 
of the airfoil, parameters coi*responding to laminar flow were used" 
from the leading edge to the point of transition,, and parameters 
corresponding to turbulent flow were used in the equations from the 
point of transition back to the trailing edge. 

Transition from laminar to turbulent flow is dependent on such 
factors as Mach number, Reynolds number, pressure gradient, stream 
turbulence, and surface roughness. Determination of transition for 
a given body and flow conditions is therefore difficult. Lees (refer- 
ence 6 ) and Schlichting (reference 7 ) have investigated the effect of 
Mach number , and velocity gradient, respectively, on the stability of 
the laminar boundaiy layer by assuming vanishingly small disturbances . 
In the absence of a stability theory which accounts for. both velocity 
gradients and compressibility, it was necessary to combine the work 
of these references. The boundary -layer thickness for neutral 
stabilily was considered to be the value for a flat plate in incom- 
pressible flow multiplied by factors to correct for Mach nuniber and 
pressure gradient. The criterion used to estimate transition in the 
present investigatipn was consideration of the n©utral stability of 
the L^nar boundary layer. The notation N = 1 is used to denote 
transition occurring when the boundary -layer momentum thickness 
reaches the value for neutral stability. If the airfoil- is well 
faired and in a stream of low initial turbulence, transition need not 
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occiir vintil the hoTjndary layer Is thicker. For exangDle^ the transition 
ctirve W = 2 corresponds to transition at a point where the hoimdary- 
layer momentiua thickness is twice that for neutral stability. 

Before a more extended treatment is possible, the stability 
theory must be developed to include the combined effects of compressi- 
bility and pressure gradient. A more complete analysis will inquire 
an investigation to determine the existence of transition regions, or 
shock waves (and hence changes in the pressure distribution) due to 
transition or separation. In the. pi*esent paper values of the momentum 
thickness of the laminar boundary layer wej:e- calci0.ated along the 
surface of, the airfoil and compared with the corresponding values of 
momentum thiclsness for nexiti’a! stability. The intersection of the 
curves through these points gives a possible point of instability 
for a given Mach number aiid Beynolds number. ^ ' The stability crlterions 
indicate that the critical part of "^e boxaidery layer from consideration 
of possible transition is well forward on the cm’ved airfoil. 

The moraent.mn thickness increases in passing through the trailing 
shock. Approximate corrections for this effect have been supplied to 
the authors by ^^r. Weal Tetervin of the Physical Eesearch Division. 

The assumptions are that the boundary -layer momentum equation applies 
and that the lengtli over which the pressure rise takes place oh the 
surface is so short that the skin friction can be neglected. The 
correction becomes: 



where 6 is the momentum thictaiess, V is the velocity tangential 
to the boundary layer, H is the average value for compressible 

‘^av 

flow of the ratio of displacement thickness to momentum thickness 
across the shock, and p is the density .at the edge of the. bovuidary 
layer. The subscripts 1 and 2 refer to conditions before and after 
the shock, respectively . 

In the actual case, the interaction of the shock wave, and 
boundary layer resvilts in an increase iii the momentum thickness, an 
increase in the displacement thickness, and a change in the velocity 
profile and may cause flow separation. The effect of separation is 
to reduce the pressure drag and increase the viscous drag. At low 
Reynolds nmbers there may be considerable separation resulting in a 
total drag less than the theoretical shock drag (reference 8) . The 
method used was chosen for lack of a' more exact method. ' 
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PRESEWTATIOW OF FIGURES 


Figuare 1(a) gives the viscous -drag coefficient of a flat plate 
as a function of Reynolds nuniber based on chord. The transition 
c\jrves for a Mach nuniber of 1*35 6 .re nunibered according to the 

ratio of "boundary -layer mcanentua tliiclsness at transition to the 
bounc ary -layer thickness for neutral stability (W = 2 to 10 ). 

Figures 1(b) and 1(c) give the corresponding curves of viscous dreg 
at Mq = 1-35 a 6 -percent -thick circular-arc airfoil and 6 -percent- 
thick double -vedge aiif oil, respectively. The transition occurs at 
a very much higher Reynolds number for the curved airfoil than for 
the flat plate or wedge, airfoil. As liie Reynolds nximber is increased, 
an abDnipt rise in the drag for the circular-arc airfoil is noticed. 

The reason for "the sudden drag rise is ihat the ratio of local 
boundary -layer thickness to the thickness for local neutral stability 
(as given by Schlichtlng 's tlieory) reaches a maxiToum well forward on 
this, airfoil BO tha,t the critical part of the airfoil (as regards 
stability or transition) is also well forward. ,A complete 'vring would 
not necessarily experience a sharp .drag rise as the Reynolds nuniber 
is increased sixice transition may occur at different Reynolds numbers 
over different sections. Parts of the transition cui-ves for the 
wedge are dashed to indicate possible theoretical . errors introduced 
by the sudden expansion at the rcidchord. The actual curves may be 
somewhat to the left of the dashed curves (closer to flat -plate 
conditions). 

Figures 2(a) and 2(b) give the viscous -drag coefficient at 
Mq = 1.6 for the 6 -percent- thick circular-arc and wedge airfoils, 

respectively. These curves show again that transition occurs at a 
much higher Reynolds number on the circular arc than on the wedgej 
however, transition on both types of airfoil occurs at a lower 
Reyn'olds nvmiber at I'Iq = 1.6 than at = I. 35 . (See fig. 1.) 

Figures 3 (S’) and 3("b) compare the viscou.s-drag coefficients of 
the airfoils at = 1.35 an^ a"?" ^ respectively, for a 

constant transition number W = 5* These figures indicate that the 
circular -arc airfoil may have a much lower viscous -drag coefficient 
than ihe wedge over a certain range of Reinolds numbers . 

The pressure -drag coefficients are determined to be as follows 
(from reference 4 )'. 
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Airfoil 

(6 percent thick) 

PresBure*drag coefficient! 

=. 1.35 

= IT6O 

Circuleir arc 
Double wedge 

0 .0218 
.0160 

0 .0158 
.0116 


IThen these pressure -dreig coefficients are added to the vi8Co^xo- 
drag coefficients of figures 3(a) and 3(^) (or in general to figs. 1 
and 2), curves are obtained for the total -drag coefficients. See 
figvires 4(a) and 4Cb). For W - 5 the circular-arc airfoil has at 
least sli^tly more drag than tlie \TOdge for all Reynolds numbers 
investigated at Mq =1.35; however, the analysis shows that the 

wedge may have the more drag over a certain Reynolds number range 
at Mq = 1.6. The actual comparison obtained depends on the value 

of. N used for the analysis. Since . N is a function of many 
parameters, such as surface finish and stream turbulence, it is 
difficiU-t to assigQ, N a proper value for a given profile. Some 
information concerning the effect of stream turbulence on transition 
at low speeds is given in reference 9 • 

This paper demonstrates the need for including both Mach number 
and Reynolds nuriber, as well as subh factors as stream turbulence 
and surface finish where possible, in papers of experimental work 
since these factors may influence the interpretation of the data for 
f^^ll-scale application. 

.It must be kept in mind that the present paper has compared the 
drags of two airfoil shapes of the same thickness ratio . A comparison 
of airfoils giving the same structural strength or stiffness wotild be 
more favorable to the curved airfoil. 


CONCLUDING REMARKS 


An analysis has been presented which semres as a preliminary 
stuciy of -the total (viscous and pressure) drag of supersonic airfoil 
sections at zero lift. .Within the limitations of the present paper 
certain conclusions have been drawn: 

The relative presstire drags of airfoils at supersonic speeds is 
in general different from the relative total drags even at zero lift. 
The airfoil shape for minimum drag varies wi-Ui Reynolds number, Mach 




KACA TN No. I 37 I 


7 


nmter, ttirljiaehce, surface finish, and other factors, and is not 
necessarily the shape that wovild give minimum theoretical pressxire 
drag. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley 'Field, Va., May 9 , I 9 U 7 
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